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The reesults of fits to RHIC-130 results, including multistrange hadrons, are shown to agree 
only with the model of an exploding QGP fireball. 
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1. Historical background 

The topic we discuss today, production of hadrons in statistical hadronization in high 
energy heavy ion collisions, has been the subject of several diploma and doctor thesis 
at the University Frankfurt in late 70’s and early 80’s. At the time, in a field void of 
any experimental result, I competed with a Hungarian hadrochemistry group lead by 
the chair of this discussion session, Prof. J. Zimanyi. 

We both passed through natural evolution stages. We were at first considering 
the equilibrium particle abundances expected to arise when Mr. EquilibriX holds in 
his hand hot and dense hadron matter fireball, which evaporated these particles. This 
work was followed by the development of kinetic theory of strangeness production, 
which was precipitated by the finding of the Budapest group, that light quark 
interactions were not fast enough to produce strangeness abundantly. 

We learned to appreciate that the yields of newly produced quarks were not 
necessarily given by the chemical equilibrium statistical model. The physics reason 
which lead us to explore the kinetic theory was relatively short available reaction 
time in relativistic heavy ion collisions, in general not allowing for the achievement of 
an equilibrium particle abundance. These remarks apply to the case that the dense 
matter fireball is made of hadrons, or deconfined quarks. Our work was extended to 
compare these two so different matter phases. 

This kinetic theory work has inspired the chemical non-equilibrium study of 
hadronization of quark-gluon plasma (QGP), and we were able to show that in 
a sudden hadronization process specific properties of QGP will be imparted in 
hadronization on particle yields. All these developments were then summarized in 
a relatively widely known theoretical review |jj . 

When the first multistrange hadron experimental results arrived 10 years ago, the 
analysis applied introduced the chemical non-equilibrium Q , which then was generally 
accepted as being an important element of data analysis. This analysis addressed solely 
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the strange particle abundances. When results for non-strange hadrons were included 
and were analyzed, it became necessary to expand the chemical-non-equilibrium 
approach, allowing the abundances of light quarks to vary 1 - 0 - In the analysis 
of the 158A GeV Pb collisions with Pb targets, we discovered that the light quark 
overabundance converged towards the maximum value allowed, at which point the 
entropy content of the final hadronic state is maximized [^| . 

In short, we expected chemical nonequilibrium, excess of entropy (hadron 
multiplicity), and abundant strange antibaryons. Experiments, too numerous to cite, 
but present in this volume, found these interesting signatures of new physics, just 
as predicted. So why are we having this discussion? As we evolved in our detailed 
understanding of the subject, we could see in our research field the arrival of new 
people who were often fast, but sometimes slow, in repeating our evolutionary steps. 

This is not the place to give those who followed the credit for not reading our 
work. However, we need to deal with the current situation that our early rather trivial 
work of 1979-81 is today taken seriously again, and is used as basis of data analysis and 
interpretation. Even if this happens without citation, this is not changing the problem 
we are facing: if this work is now widely reported, is this ‘majority’ vote suggesting that 
the very extensive theoretical development carried out by many prominent theorists 
over the past 20 years has gone in the wrong direction? Clearly, this is not the case, 
so we have at least to make an effort to end this rather absurd situation. 

The principal speaker (K. Redlich) has presented his views, which we address in 
technical manner in our main contribution to this volume This contribution is a 
qualitative discussion of approaches and results in the study of the experimental data. 

2. Strangeness Equilibration 

Collision of small elementary systems produce colored partons linked by strings, which 
‘snap’ giving birth to particle pairs, and strangeness is in p-p interactions produced 1/5 
as often compared to the light flavor u or d. The reasons this happens are understood 
in the realm of string tension and particle tunneling. Such results form the basis of 
simulations of nuclear collisions with independent string mechanics. 

The formation of a deconfined region of space-time is at origin of another effective 
mechanism of strangeness production. Beside (confined or deconfined) quarks, we now 
also have gluons in the hot gas, and glue-glue collisions are believed to be the source 
of the abundant strangeness observed in the experiment. The string mechanism uses 
the chromo-electric field between color charges, the gluon is the transverse quantum 
of the chromo-fields. Thus, though in both cases mechanisms specific to the glue filed 
are invoked, both mechanisms are complementary. 

It is a mere coincidence that given the strength of the string tension (1 GeV/fm), 
the mass of the strange quark (160 ± 50 MeV), and the initial temperatures we reach 
in the hot fireball (200-300 MeV), that in the p-p interactions and central Pb-Pb 
interactions the yield of strangeness is near to chemical equilibrium yield. A lot 
of theoretical effort went to understand this situation within kinetic models. This 
fact cannot be generalized to different collision systems without prior and thorough 
study of the kinetic production mechanism applicable. The extension of the work 
on strangeness is the study of light quark abundance. It is not guaranteed that at 
all conditions, and for all systems considered even u , d quarks are found in chemical 
equilibrium. 
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Figure 1. Strangeness chemical relaxation time r s , for a s (Mz) = 0.118 with m s (Mz) = 
90 MeV, pj°(m s ~ 200 MeV) (thick line). Hatched areas: effect of variation of strange quark 
mass by 20%. 


The deviations from chemical equilibrium can be incorporated into the momentum 
distribution of (strange) quarks and antiquarks in the deconfined matter in the form: 


/s/s 


7s" 1 A + 1 


7sA;re 


±1 -VV+mf/T 


(1) 


When 7 S —> 1 we are approaching chemical equilibrium distribution. On the right hand 
side in Eq. (jl[), we see the Boltzmann limit. This distribution maximizes the entropy 
at given particle number. When a Boltzmann momentum distribution is established, 
the phase space occupancy 7 s (t) may be still quite different from unity, as indicated it 
evolves in time t due to ongoing pair production (gluon fusion) processes, GG —> ss. 

The characteristic time required to equilibrate strangeness chemically in the 
deconfined gluon gas is shown Fig.[l]. This result suggests that depending on the 
initial temperature and subsequent evolution, the lifespan of the QGP needs to live 
at high temperature in excess of 4 fm/c in order for the strange quark flavor to reach 
full chemical equilibrium. This lifespan of the plasma phase is not generally assured, 
and thus in QGP the chemical equilibrium is not easily reached. 


3. Hadronization 

Hadrons are always the most abundant final particles observed experimentally. We 
must be able to understand and interpret the production mechanism of these particles. 
Strange hadrons are the largest particle family, and are thus particularly interesting. 
Our primary objective is to identify what changes can be expected in the yields of 
strange particles emitted from deconfined quark-gluon matter. And we must device a 
mechanism of data analysis that is sufficiently sensitive to recognize this. 

For now nearly 50 years we know that the abundances of hadrons produced in 
elementary and heavy ion collisions can be, within a factor of two, described by the 
statistical phase space. Hagedorn and before him, Fermi ||, Landau © , have 
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pioneered that idea. One can argue that when the production probability (the square 
of the quantum matrix element) of different hadronic particles is similar, they saturate 
by strength of their interactions the unitarity limit, the primary difference in particle 
yield is arising from the relative magnitude of the accessible emission phase space. 

This reduces the issue to the question, why can we use a ‘temperature’ T 
parameter to characterize the size of hadronic particle phase space? Hagedorn spoke 
of a pre-existent temperature, later he developed a bootstrap model which helped 
quantify this idea |TlJ. Another possible approach uses the quark-vacuum fluctuations 
to obtain Boltzmann-like hadron momentum distribution |12[. In nuclear collisions, 
a large number of quarks and gluons within a deconfined local domain of space time 
can approach by two body collision processes the local thermal distribution within 
a relatively short period of time. In both cases, the confining vacuum is the critical 
input factor, as it keeps together these strongly interacting quanta. In my view, the 
confining vacuum makes strong interactions different and introduces the statistical 
phase space into hadronic particle abundances. 

Two properties, both associated with gluons in the deconfined state, allow 
to remember the structure of the fireball, and eventually help us understand the 
properties of QGP. The addition of gluonic degrees of freedom is doubling the entropy 
content of quarks alone, which already benefit from the breaking of color bonds to 
increase in the deconfined state the entropy content; gluon based reactions are effective 
in producing strange quarks. The enhancement of the entropy content is seen in 
enhanced pion yields, while enhanced strangeness abundance contributes to a pattern 
of strange hadron which favors production of multistrange hadrons. 

What is the level of sensitivity in particle detection which is required for these 
observables? The entropy enhancement translates into a 30-60% increased pion yield. 
Strangeness enhancement is expected at the level of factor two-three. One would 
imagine, seen the magnitude of these effects of QGP, that we need to measure these 
quantities at the level of a few percent. Thus, we must obtain experimental results 
at this level of precision, and device an analysis method which is sensitive at this 
precision level. 

In quantitative terms, the relative number of final state hadronic particles freezing 
out from, e.g., a thermal quark-gluon source, is obtained noting that the fugacity /) 
of the i-th emitted composite hadronic particle containing ^-components is derived 
from fugacities A k and phase space occupancies 7 k' 

Ni oc e~ Ei/Ti fi = e ~ Ei/Tf JJqfcAfc. (2) 

k£i 

As seen in Eq. ®), we study particle production in terms of five statistical parameters 
T, A q , A s , 7 q , 7 S . The difference between the two types of chemical parameters, A; and 
7i, is that the phase space occupancy factor 7* regulates the number of pairs of flavor 
V, and hence applies in the same manner to particles and antiparticles, while fugacity 
A i applies only to particles, while A ” 1 is the antiparticle fugacity, as is seen in Eq. (|l|). 

To describe the shape of spectra, one needs matter flow velocity parameters, these 
become irrelevant when total particle abundances are studied, obtained integrating all 
of (the velocity deformed) phase space. In presence of strong longitudinal flow, the 
‘scaling’ in rapidity implies that when we are looking at a yield per unit of rapidity, 
we can also ignore the velocity parameter when considering particle ratios. 

The resulting yields of final state hadronic particles are most conveniently 
characterized taking the Laplace transform of the accessible phase space. This 
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approach generates a function which, in its mathematical properties, is identical to 
the hadron gas partition function: 


C 


e Ei/Tt 7iA- 

k£i 


oc In Z HG . 


( 3 ) 


This characterization of the accessible phase space does not require formation of a 
phase comprising a gas of hadrons, but is not inconsistent with such a step in evolution 
of the matter. The final particle abundances, measured in an experiment, are obtained 
after all unstable hadronic resonances ‘j’ are allowed to disintegrate, contributing to 
the yields of stable hadrons. The unnormalized particle multiplicities are obtained 
differentiating Eq. (|3]) with respect to particle fugacity. 

Are we sure that all particle yields produced predominantly in the statistical 
hadronization process? Already in 1982, we were pointing out that rarely produced 
particles such as f l, would be the first to show additional enhancement originating 
in either new mechanisms, or simply in the fact that strangeness correlation would be 
favoring strange clusters which alter statistical hadronization yields of ft, O Jisf . 

It is important that the reader understands the method how one finds such 
exceptions from the statistical hadronization pattern. Given that we have many 
available particle yields, we can look at the data fits omitting one result at a time, in 
order to see that the fit is stable, and that it does not depend on one single experimental 
result. In order to report statistical significance, this procedure must be carried out, 
it is not implemented in fit programs. 

The omission of the Q yields (and to a lesser extend Q) in the study of SPS results 
at i/snn = 17-2 GeV, had the effect of allowing for a much better convergence of the 
statistical hadronization fit. This means that within the reported f 1 (and f2) particle 
yields there are additional contributions outside of the statistical hadronization 
picture. Subsequent study of the mj_-spectra of S2 + O have in fact confirmed that 
there is an additional soft momentum component |l4|| . There are several explanations 
today of this behavior. Whatever is the true mix of production mechanisms, we simply 
do not include the fl, yields in the statistical hadronization analysis, keeping in mind 
that their large yield is a challenge for cascade model simulations. 


4. Hadron data fits 


We turn rightaway to discuss the RHIC data obtained in ^/snn = 130 GeV run. 
The experimental results we consider are the (combined where possible) data of 
STAR, PHENIX, BRAHMS, PHOBOS, for more discussion of the data origin, see 
0 - For comparison purposes, we adopted the results nearly as presented in this 
reference. However, we did not fit natural results such as 7r + /7r“ = 1. We also do 
not use the results for K*, K* since these yields depend on the degree of rescattering 
of resonance decay products. We assume, in our fits in table that the multistrange 
weak interaction cascading S —> A, in the STAR result we consider [|h|, is cut by 
vertex discrimination and we use these yields without a further correction. 

The data analysed are at the central rapidity region where, due to approximate 
longitudinal scaling, the effects of flow cancel and we can evaluate the full phase space 
yields in order to obtain particle ratios. In the last column in table |], the chemical 
equilibrium fit, the large x 2 originates in the inability to account for multistrange 
S, 5. Similar results as here shown 116 , are presented in Ref. 0 , which work 


was 


done before the multistrange hadron yields were reported. 





Strangeness, Equilibration, Hadronization 


6 


Table 1 . Fits of central rapidity hadron ratios for RHIC -^snn = 130 GeV run. Top 
section: experimental results, followed by chemical parameters, physical property of the 
phase space (energy E and entropy S per baryon number b), and the fit error. Columns: 
data, full non-equilibrium fit, nonequilibrium fit constrained by strangeness conservation 
and supersaturation of pion phase space, and in the last column, equilibrium fit 
constrained by strangeness conservation, upper index * indicates quantities fixed by 
these considerations. 



Data 

Fit 

Fit 

Fit eq 




s — s = 0 

s — s = 0 

p/p 

0.64 ±0.07 

0.637 

0.640 

0.587 

p/h" 


0.068 

0.068 

0.075 

A/A 

0.77 ±0.07 

0.719 

0.718 

0.679 

A/h" 

0.059±0.001 

0.059 

0.059 

0.059 

A/h" 

0.042±0.001 

0.042 

0.042 

0.040 

s/s 

0.83 ±0.08 

0.817 

0.813 

0.790 

S"/A 

0.195±0.015 

0.176 

0.176 

0.130 

S-/A 

0.210±0.015 

0.200 

0.200 

0.152 

K-/K+ 

0.88 ±0.05 

0.896 

0.900 

0.891 

K — / 7T~ 

0.149±0.020 

0.152 

0.152 

0.145 

Kg/h- 

0.130±0.001 

0.130 

0.130 

0.124 

n/z~ 


0.222 

0.223 

0.208 

Ti/Zr 


0.257 

0.256 

0.247 

n/ ci 


0.943 

0.934 

0.935 

T 


158± 1 

158± 1 

183± 1 

7q 


1.55±0.01 

1.58±0.08 

1* 

Aq 


1.082±0.010 

1.081±0.006 

1.097±0.006 

7s 


2.09±0.03 

2.1±0.1 

1* 

A s 


1.0097±0.015 

1.0114* 

1.011* 

E/b[ GeV] 


24.6 

24.7 

21 

s/b 


6.1 

6.2 

4.2 

S/b 


151 

152 

131 

E/S[MeV] 


163 

163 

159 

X^/dof 


2.95/(10-5) 

2.96/(10-4) 

73/(10-2) 


The chemical equilibrium fit yields E/S = 159 MeV < T = 183 MeV, and thus 
in order to satisfy the constraints arising from the Gibbs-Duham relation: 


E 

Is 


PV 

~s~ 


T + ST , 



Mb = 3T In A q , 


(4) 


the emission volume must be unrealistically large, since the pressure is much smaller 
than the energy density E/V. This inconsistency is to best of my knowledge well 
known to the promoters of the equilibrium hadronization model, yet they do nothing 
to resolve it. 

On the other hand, in the chemical nonequilibrium fits with 7 s ,7 q > 1 there 
is good agreement with experimental data. Results are as would be expected for a 
rapidly expanding fireball, they converge to the maximum allowed value: 


7q = 7q = e 


m 7r /2Tf 


(5) 
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which corresponds to the maximum entropy in the pion gas. The large values of 
7 q > 1 is consistent with the need to hadronize the excess entropy of the QGP. The 
volume problem we addressed above is absent since E/S = 163 > T = 158 MeV. The 
value of the hadronization temperature T = 158 MeV is below the central expected 
equilibrium phase transition temperature, and this hadronization temperatures at 
RHIC is consistent with sudden breakup of a supercooled QGP fireball. We conclude 
that as has been predicted 15 years ago the inclusion of the yields of multistrange 
antibaryons in the RHIC data analysis, along with allowance for chemical non¬ 
equilibrium (7 / 1 ), has the sensitivity to discriminate the different reaction scenarios. 

The value of the thermal energy content E/b = 25 GeV, seen in table [|, is in 
very good agreement with expectations once we allow for the kinetic energy content 
associated with longitudinal and transverse motion. The energy of each particle is 
‘boosted’ with the factor 7 " coshj/||. For u_l = c/\/3, we have 7 " = 1.22. The 
longitudinal flow range is about ±2.3 rapidity units, according to PHOBOS results. 
To obtain the energy increase due to longitudinal flow, we have to multiply by the 
average, / dy\\ cosh y\\/y\\ —► sinh(2.3)/2.3 = 2.15, for a total average increase in energy 
by factor 2.62, which takes the full energy content to E v /b ~ 65 GeV as expected. 
This consistency reassures that we have a physically relevant fit of the data, which 
respects overall energy conservation. 

The strangeness content, s/b = 6 , is slightly below the expected equilibrated 
QGP phase space, which would have yielded 8.6 strange quark pairs per baryon at 
A q = 1.08. Thus 7 ^ gp = 6 / 8.6 ~ 0.7, which is greater than the value 0.5 a similar 
discussion yields for the SPS energy range. Using the fitted value 7 PG = 2.1, we find 
at RHIC like at SPS, 7 PG / 7 ? GP — 3, which is the correct phase space size ratio for 
HG and QGP at these conditions. 

5. Final remarks 

There is ample evidence, looking at particle spectra at RHIC, that the evolution of 
the system is very rapid. We have pioneered this reaction model in the analysis of the 
SPS data. We have in particular pointed out that the explosive disintegration of the 
QGP accompanied by supercooling below equilibrium phase transition temperature 
is a natural phenomenon flg|| . Thus, in the hadronization model, we study this is 
the natural background against which we make decisions how to proceed and how to 
interpret the results obtained in analysis of particle spectra. 

In sudden hadronization, V HG ~ yQG p ) the growth of volume is negligible, 
T QGP ~ T hg , the temperature is maintained across the hadronization front, and 
the chemical occupancy factors in both states of matter accommodate the different 
magnitude of the particle phase space. In this case, the QGP strangeness when 
‘squeezed’ into the smaller HG phase space results in 7 PG / 7 ^ GP — 3, which is of 
the same magnitude as the unfrozen color degeneracy. This theoretical expectation is 
observed both at the top SPS and at the RHIC-130 run. 

For the top SPS energy range, this interesting result could be ignored, since 
accidentally 7 PG ~ 7 PG - For this reason hadronization at SPS energy can be also 
modeled in terms of an equilibrium hadronization model: The pion enhancement 
associated with the high entropy phase can be accommodated by use of two 
temperatures, one for the determination of absolute particle yields, and another 
for determination of the spectral shape. Such an approach has similar number of 
parameters, as the method presented here, the only inconsistency (with HBT) is the 
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However, already at RHIC-130 the hadron phase space occupancy for strangeness 
is significantly larger than for light quarks, see table [j]. It is the inclusion of the yields 
of multistrange antibaryons in the RHIC data analysis, which leads to clear resolution 
of the hadronization mechanism. The specific per baryon strangeness yield at RHIC 
is an order of magnitude greater than at SPS. 

We see at RHIC considerable convergence of the hadron production around 
properties of suddenly hadronizing entropy and strangeness rich QGP. There 
is overwhelming theoretical and experimental evidence that the chemical non¬ 
equilibrium yields of hadrons are required to understand with the necessary precision 
the nature of the quark-gluon fireball formed in these interactions. 
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